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"> ; ABSTRACT 

Results are presented on the X-ray properties of 9 high-redshift (1.2 < z < 3.4) 
(Nj , radio-loud quasars (RLQs) observed by ASCA (10 observations) and ROSAT (11 

observations, for a subset of 6 quasars). New ASCA observations of S5 0014+81 (z 
= 3.38) and S5 0836+71 (z = 2.17) and ROSAT observations of PKS 2126-158 for 
ON ■ which results were never presented elsewhere, are included. 

A simple model consisting of a power law plus cold, uniform absorption gives ac- 
ceptable fits to the spectra of all sources. The ASCA spectra of the 6 brightest objects 
show evidence for absorption in excess of the Galactic value at a ^ 99% confidence 
level. Comparison with the ROSAT data, suggests that absorption has significantly 
varied (AA H ~ 8 x 10 20 cm -2 ) in the case of S5 0836+71, on a time-scale of ~ 0.8 



yr in the quasar frame. For the remaining 5 sources for which ROSAT spectra were 
available, the two instruments gave consistent results and the data were combined 
yielding unprecedent spectral coverage (typically ~ 0.4-40 keV in the quasar frame) 
for high-z quasars. This allows to put severe limits on several different descriptions 
of the continuum (e.g. broken power law, bremsstrahlung, reflection component). No 
Fe Ka emission line is detected in any of the ASCA spectra. An absorption edge 
consistent with Fe Ka at the quasar redshift is marginally detected in S5 0014+81. 
Possible origins for the observed low energy absorption are discussed. In particular, 
contributions from the molecular clouds and dust present in our Galaxy (usually dis- 
regarded) are carefully considered. In the light of the new results for S5 0836+71 and 
S5 0014+81, absorption intrinsic to the quasars is considered and discussed. 

The average slope obtained from the 8 ASCA spectra in the observed ~ 0.5-10 
keV energy band is < To.s-iokcV > — 1.61 ± 0.04, with a dispersion o"o.5_iokcv — 0.10 
± 0.03. The average photon index in the observed 2-10 keV band, where the effect of 
absorption is negligible, is < r2-iokcV > — 1-53 ± 0.05, with a dispersion o"2_iokev ~ 
0.12. Furthermore, the implications of the present results on the calculations of the 
contribution of quasars to the cosmic X- and 7-ray backgrounds (XRB and GRB) are 
briefly discussed. 

Subject headings: galaxies: active — quasars: general — X-rays: galaxies 
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1. Introduction 

Quasars are the most powerful objects in the 
whole Universe. This is especially true in the X- 
ray band, where luminosities can reach ~ io 47 ~ 48 
erg s _1 . However how quasars produce such a large 
amount of energy remains a challenging astrophys- 
ical problem. Certainly, because of their extreme 
conditions, quasars provide a powerful test for mod- 
els of emission mechanisms of active galactic nuclei 
(AGNs) (Rees 1984). Quasars show strong contin- 
uum emission over the entire electro-magnetic spec- 
trum, from radio through the X- and even 7-ray 
region (Sanders et al. 1989, Elvis et al. 1994a, 
Thompson et al. 1995). Optically selected samples 
of quasars indicate that ~ 90% of them are radio 
quiet (RQQs) and ~ 10% arc radio loud (RLQs). 

X-ray quasar spectral observations are crucial for 
two main reasons: X-rays carry a large amount of 
the total quasar luminosity. Secondly, as demon- 
strated by observations of rapid X-ray variability, 
X-rays originate from the innermost regions of the 
quasar (Mushotzky, Done & Pounds 1993). Most X- 
ray spectral observations have included mainly low 
redshift (z < 1) quasars; the poor energy resolu- 
tion generally limited the analysis to a simple pa- 
rameterization of the spectrum with a single power 
law. In the ~ 0.1-4 keV energy range, previous 
Einstein IPC and ROSAT PSPC observations have 
shown that RLQs have significantly flatter X-ray 
spectra than RQQs (Wilkes & Elvis 1987, Brun- 
ner et al. 1992) and that, for a given optical lumi- 
nosity, RLQs are on average ~ 3 times brighter in 
X-rays than RQQs (Zamorani et al. 1981). Stud- 
ies at higher energies (~ 2-10 keV) with EXOSAT 
and Ginga have confirmed the dichotomy, with a 
clear correlation between spectral index and radio- 
loudncss (Williams et al. 1992, Lawson et al. 1992). 
Whether differences in the observed X-ray proper- 
ties should be attributed either to intrinsically dif- 
ferent properties of the sources or inclination effects 
and/or host galaxy properties is not yet well under- 
stood. Also selection effects and/or complex spec- 
tral structures (e.g. soft-excess emission, ionized ab- 
sorption) may complicate the correct interpretation 
of the data (Halpern 1984, Comastri et al. 1992, 
Fiore et al. 1993). 

At high rcdshifts (^1), spectral information is al- 
most absent for RQQs, and scarce for RLQs, since 
soft X-ray observations (mainly from the ROSAT 
PSPC) have allowed only a poor determination of 
the spectral slopes for only a small number of ob- 
jects, mostly RLQs (Bechtold et al. 1994, Elvis 
et al. 1994b, hereafter E94). It is not yet clear 
whether quasars do exhibit spectral evolution. This 
is a fundamental question which has a direct impact 



on quasar formation models. A remarkable result 
has been the discovery made with the PSPC that 
at least some of the high-z RLQs have low energy 
cut-offs possibly due to absorption along the line of 
sight (E94). Preliminary results of ASCA observa- 
tions of high-z RLQs have indeed already confirmed 
the low energy cut-off in two of them (Serlemitsos 
et al. 1994) and discovered a third probable case 
(Siebert et al. 1996, hereafter S96). Comparison 
with two high S/N spectra of two RQQs at z ~ 
1 (Nandra et al. 1995) indicates that RQQs are 
steeper than RLQ even at z > 1. However, the 
number of quasars observed so far is too small to 
draw any reliable conclusion. 

This paper presents a comprehensive and uniform 
study of ASCA observations of a sample of 9 RLQs 
with 1.2 < z < 3.4. Whenever possible, the ASCA 
results are compared and combined to ROSAT spec- 
tra extracted from the public archive. Extensive 
search for Fe K emission lines, high-energy excesses 
("hard tails") and alternative models are presented. 
The possible origin of the apparently common ex- 
cess absorption found in the data is discussed in 
the light of two newly discovered RLQs with such 
feature. Finally, the impact of these new measure- 
ments on the cosmic high energy background radi- 
ation is briefly discussed. 

In the following, Ho = 50 km s _1 Mpc -1 and 
qo=0 are assumed throughout. 

2. Observations and Data Analysis 
2.1. The Sample 

The quasar sample consists of all objects (9) ei- 
ther from PI or archival ASCA observations avail- 
able before 1996, January 1 st . The purpose was 
to analyze data of a reasonably large number of 
quasars to be able to address for the first time sta- 
tistically the X-ray properties of the class. 

A total of 10 ASCA pointed observations were 
collected, with S5 0014+81 observed twice. For 
five quasars, 11 ROSAT PSPC observations were 
retrieved from the archive and the source spectra 
re-analyzed to ensure a uniform and consistent anal- 
ysis within the sample. A considerable part of our 
analysis reproduces in part previous work on indi- 
vidual sources (see §4). The present analysis, how- 
ever, differs from single object studies as: a) it pro- 
vides a uniform analysis of the quasar sample, b) it 
makes use of the most recent calibrations (partic- 
ularly important for those observations performed 
during the Performance and Verification phase of 
ASCA) and c) it compares on a uniform basis ASCA 
and ROSAT spectral results. It is worth pointing 
out that ASCA results on S5 0014+81 (z=3.38) and 
S5 0836+71 (z=2.17) are new. The relevant data for 
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the whole sample are given in Table 1. 

2.2. ROSAT Data Reduction 

A sub-sample of six of the 9 quasars were ob- 
served on-axis with the Position Sensitive Propor- 
tional Counter (PSPC) (Pfcfferman et al. 1987) on 
board the R OS A T Observatory (Triimper 1983) be- 
tween 1991 and 1993. The relevant data for the ob- 
servations are listed in Table 2. As indicated in the 
Table, two new PSPC observations of PKS 2126-158 
are presented, which almost double the total num- 
ber of counts available for that source. The PSPC 
has an energy bandpass in the range 0.1-2.4 keV 
with an energy resolution of AE/E <~ 0.5 keV at 1 
keV. Source spectra were extracted from circles of <~ 
90" -200" centered on the sources, and background 
spectra were taken from annuli centered on the 
sources or from circular regions uncontaminated by 
nearby sources. Source and background counts were 
corrected for telescope vignetting. Data preparation 
and analysis were performed using the JAN95 ver- 
sion of the EXSAS/MIDAS software package (Zim- 
merman et al. 1993). Spectral analysis was per- 
formed using the version 8.50 of the XSPEC pro- 
gram (Arnaud, Haberl & Tennant 1991). 

2.3. ASCA Data Reduction 

ASCA has two gas imaging spectrometers (CIS) 
and two solid-state imaging spectrometers (SIS) 
(Tanaka, Inoue & Holt 1994). The energy resolution 
of the CIS and SIS are AE/E - 0.15 keV and - 0.05 
keV at 1 keV, respectively, which is about 3 and 
10 times better than the ROSAT PSPC. The SIS 
was operated in 1-CCD, 2-CCD or 4-CCD modes, 
depending on the observation (see Table 3). Only 
chip n.l of SIS0 and chip n.3 of SIS1 were used in the 
analysis of the SIS data, except for the AOl obser- 
vation of S5 0014+81 where also chip n.2 of SIS1 was 
used because the source photons were spread equally 
over the two detectors. Following a software-related 
problem on board ASCA, the data collected from 
the GIS3 during the observation of PKS 0537-286 
were damaged. They could not be recovered and 
therefore were excluded from the analysis. All ob- 
servations were performed in FAINT mode and were 
corrected for dark frame error (DFE) and echo un- 
certainties (Otani & Dotani, 1994). The data were 
selected according to standard (rather conservative) 
criteria, i.e., when the angle between the target and 
the Earth's limb was > 5° (SIS and CIS), when the 
geo-magnetic rigidity was > 7 GeV/c (GIS) and > 
8 GeV/c (SIS) and when the angle between the tar- 
get and the day-night terminator was > 25° (SIS). 
Rise-time rejection was applied to the GIS data, and 
hot and flickering pixels were removed from the SIS 
data. Telemetry drop-outs and spikes were excluded 



from the light curves for each instrument. 

Source counts were extracted from circles cen- 
tered on the sources of 6' for the GIS and 3' for the 
SIS. For the SIS instruments, the background spec- 
tra were obtained from the edges of the same CCD 
chip. The use of SIS blank sky files for the back- 
ground yielded spectral results within the errors re- 
ported in the following analysis. For the GIS instru- 
ments, we found that the two standard background 
substraction methods (backgrounds extracted from 
the blank sky files or locally) gave systematically 
different results for the weakest sources. There- 
fore, a non-standard background extracted from the 
blank sky files was adopted. Detailed explanation 
and justification for this choice is given in the Ap- 
pendix A. 

The relevant data for the ASCA observations are 
given in Tabic 3. Data preparation and spectral 
analysis were performed using version l.Oh of the 
XSELECT package and version 8.50 of the XSPEC 
program (Arnaud et al. 1991). 

3. Results 

3.1. ROSAT Temporal and Spectral Analysis 

The data were first binned in 400 s time intervals, 
as suggested by the wobble period of the telescope, 
and light curves were plotted to evidence variabil- 
ity. However, no significant variation was detected. 
Source spectra were thus accumulated for each ob- 
servation, And binned with a signal to noise ratio 
from 4 to 13, depending on the source statistics. 
Three quasars (PKS 0438-436, S5 0836+71 and PKS 
2126-158) had multiple observations. Comparisons 
between different observations indicate a clear flux 
variation for S5 0836+71, and possibly for the other 
two quasars. These are discussed in detail in §4 
where results for individual sources are presented. 
After preliminary analysis, no significant spectral 
variability was detected between the multiple ob- 
servations. These were, therefore, fitted simultane- 
ously, tying together the fit parameters, but with 
the normalizations free to vary. 

Fits were performed using a single absorbed power 
law model with absorption cross-sections and abun- 
dances from Morrison & McCammon (1983). The 
resulting spectral photon indices T, column densities 
Nn and normalizations are given in Table 4. The 
two-dimensional \ 2 contour plots in the parameter 
space Afj-r are shown in Fig. 1, together with the 
ASCA contours (see §3.2). The present re-analysis 
yielded results consistent with previous measure- 
ments (see references in Table 2). It is worth point- 
ing out that the absorption in excess of the Galac- 
tic value previously reported for PKS 0438-436 and 
PKS 2126-158 (Wilkes et al. 1992, E94) is con- 
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firmed by the present analysis (see Fig. 1 or Table 
4), which makes use of two new PSPC observations 
of PKS 2126-158 almost doubling the total counts 
from that source. The weighted mean photon index 
for the PSPC sample is < ro.i-2.4kcV > — 1-53 ± 
0.06. 

More complex models were not attempted since 
a) the main purpose of the present analysis of 
ROSAT observations is to compare and, whenever 
possible, combine these data with the ASCA data 
and b) in all cases, a single absorbed power law 
model provides an acceptable description of the 
spectra. More complex models might, however, be 
found in the references given in Table 2. 

3.2. ASCA Temporal and Spectral Analysis 

Source plus background light-curves were accu- 
mulated for each source, and none of these indi- 
cated significant variability. This is not surprising 
given the large statistical scatter of the data due to 
the low counting rates. GIS and SIS spectra were 
binned with more than 20 counts/bin between ~ 
0.7-10 keV and - 0.5-10 keV, respectively. The 
matrices used were the gis[23]v4-0.rmf released in 
June 1995 for the GIS and the "rspl.lalphaPI" ma- 
trices released in October 1994 for the SIS. Since 
the spectral parameters obtained by separately fit- 
ting the four detectors were all consistent at a ~ 
90% confidence level, the data were fitted simul- 
taneously from all four instruments with the same 
model, tying the fit parameters together but allow- 
ing the relative normalizations of the four datasets 
to vary. 

The spectra were first fitted using a single ab- 
sorbed power law model with all parameters free to 
vary. The resulting best-fit parameters are given in 
Table 5, together with the absorbed 2-10 keV flux 
and intrinsic 2-10 keV luminosity derived from the 
fits. In each case, a single absorbed power law model 
provides an acceptable description of the spectra. 
The two-dimensional x 2 contour plots in the pa- 
rameter space -/Vh-T are shown in Fig. 1, together 
with the ROSAT contours. Contours representing 
the 68%, 90% and 99% confidence limits for two 
interesting parameters are indicated for the simul- 
taneous fit of the GIS and SIS data. The 90% confi- 
dence contours obtained from separately fitting GIS 
and SIS are indicated as well. They clearly show 
that the spectral parameters obtained from the GIS 
always agree with the SIS results at least at a 90% 
confidence level. The elongated shape of the GIS 
contours in the direction of low column densities is 
a consequence of the reduction of the GIS effective 
area at low (E < 0.8 keV) energies, where the ef- 
fect of the absorption is highest; to be compared 
with the SIS detectors sensitive down to E ~ 0.4- 



0.5 keV. The much smaller phase space occupied 
by the ASCA x 2 contours compared with ROSAT, 
in particular for S5 0014+81 and NRAO 140, indi- 
cate the higher capabilities of ASCA to constrain 
both the photon index (because of the larger energy 
band) and column density (because of the higher 
statistics). 

The contours also clearly show that at least 6 
quasars: S5 0014+81, NRAO 140, PKS 0438-436, 
S5 0836+71, PKS 2126-158 and PKS 2149-306, have 
Nn values larger than the Galactic one, at a more 
than 99% confidence level. For the three Parkes 
sources and NRAO 140, the present results confirm 
previous findings obtained with the same ASCA 
data (Serlcmitsos et al. (1994), S96, Turner et al. 
1995). The most important results are instead for 
S5 0014+81 and S5 0836+71 for which excess ab- 
sorption is detected for the first time, at a high (> 
99.99%) confidence level (Fig. 1). The previous 
analysis of the ASCA observation of S5 0014+81 by 
Elvis et al. (1994c) didn't show any excess absorp- 
tion as only the GIS data had been analyzed, thus 
sensitivity at low energies to detect the excess ab- 
sorption was lost. A striking result from the present 
analysis comes from comparing the ROSAT and 
ASCA contours for S5 0836+71 (Comastri, Cappi & 
Matsuoka 1996) which indicates a large variation in 
the absorption column density of ~ 8 x 10 20 cm~ 2 
(see §4.6 for details). 

It should be noted, however, that although the 
X 2 contour plots of Fig. 1 show that the ASCA and 
ROSAT parameters are always consistent within 
their 90% confidence levels (except for S5 0836+71), 
the same contours also indicate that the best-fit col- 
umn densities obtained with ASCA are systemat- 
ically larger than the ROSAT values. There are 
several possibilities to explain this apparent contra- 
diction: (1) all quasars experienced a real increase 
in the absorption column densities (between ~ 6— 
15 x 10 20 cm -2 ); (2) all the excess absorption is 
due to SIS calibration uncertainties at low energies; 
(3) the absorption measured by ASCA (corrected 
for the calibration uncertainties estimated below) 
is real and the difference between the ASCA and 
ROSAT results is due to a) the increased statis- 
tic obtained with ASCA or b) the fact that the 
absorption is more complex than the adopted one 
(i.e. a warm absorber). Absorption variability (1) 
can be disregarded as there is no reason to expect 
such increase in all quasars. In order to quantify 
any systematic instrumental effect (2) of the SIS at 
low energies, a series of tests described in detail in 
Appendix B have been applied. Results from this 
study indicate that part (~ 2-3 x 10 20 cm -2 ) of 
the excess absorption column density measured in 
the quasars can indeed be attributed to remaining 
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SIS calibration uncertainties. However it is very un- 
likely that the effect is all instrumental since even 
after considering a conservative systematic error of 
~ 3-4 x 10 20 cm~ 2 (estimated by us in the Ap- 
pendix B and, independently by Hayashida et al. 
1995), the measured column densities are still sig- 
nificantly higher than the Galactic values. Also in 
the light of the fairly good agreement between the 
ROSAT and ASCA results, in particular for PKS 
0438-436 and PKS 2126-158 (see Fig. 1), we arc in- 
clined to interpret these excess absorptions as real. 
Either the improved ASCA signal to noise ratio or a 
complex absorber, therefore, are the most plausible 
and will be further considered below (see §3.3.2). 

The weighted mean photon index obtained for 
the whole sample in the observed ~ 0.5-10 keV is 

< r .5-iokcV > — 1-59 ± 0.01 with a dispersion 
00.5-iokeV — 0.10 when Nh is left free to vary 
and 1.45 ± 0.01 with a dispersion of 0.08 when 
Nh is fixed at the Galactic value. However, fits 
with the absorption fixed at the Galactic value (Ta- 
ble 5) are not acceptable for these quasars because 
clear depressions appear at low (< 1 keV) energies 
in the residuals of SIS and, in some case, of GIS. 
The mean slope and dispersion have also been com- 
puted, jointly, using a maximum-likelihood tech- 
nique (see Maccacaro et al. 1988) which has the 
advantage, over a simple algebraic mean, of weight- 
ing the individual photon indices according to their 
measured errors assuming a Gaussian intrinsic dis- 
tribution of spectral indices. The results obtained 
in this way are < r .5-iokcV >— 1-61 ± 0.04 and 
co.5-iokeV = ®AQ ± 0.03, where the confidence in- 
tervals are the 68% level for two interesting pa- 
rameters. Fitting the data only between 2-10 keV 
where the effect of the absorption is negligible gives 

< r 2 -iokcV > ^ 1-53 ± 0.05 and <7 2 -iokcV ^ 0.12. 
The ASCA spectra were also investigated for the 

presence of Fe Ka emission lines. Given the high 
rcdshifts of the quasars, Fc Ka lines are expected 
at energies between ~ 1.5 and 3 keV, where GIS 
and SIS effective areas and resolution are highest. 
No significant Fe Ka emission lines were detected, 
with upper limits for the equivalent width of a nar- 
row Gaussian line at 6.4 keV ranging between 40 
and 415 eV (quasar frame, Table 5). Very similar 
upper limits were obtained for a line emitted at 6.7 
keV in the quasar frame. It is emphasized, how- 
ever, that since no constraint can be given on the 
lines widths, the upper limits strongly depend on 
the assumed widths of the lines and become a factor 
two worse than given here if the widths are as large 
as 0.5 keV. Next, spectra were searched for neu- 
tral Fe K edges at 7.1 keV (quasar frame) and, for 
ionized Fe K edges, at the mean energy of 7.8 keV 
observed in Seyfert 1 galaxies (Nandra & Pounds 



1994). Whenever these spectral features were ex- 
pected in the energy range of 1.8-2.4 keV, a conser- 
vative systematic error of 10 eV (for the lines) and 
At = 0.05 (for the edges) were added to take into 
account the uncertainties of the ASCA response at 
these energies (Table 5). No absorption edges were 
detected except for S5 0014+81. For this quasar, 
the inclusion of an absorption edge at E ~ 1.62 ± 
0.1 keV (E ~ 7.1 ± 0.4 keV in the quasar frame) 
with depth r ~ O.^t^;} 2 corresponds to a A% 2 
~ 9 which is significant at more than 90% confi- 
dence. The detection is also supported by the fact 
that the edge seems to be present in both obser- 
vations (October 93 and October 94), at 68% and 
99% confidence level, respectively. The equivalent 
hydrogen column density derived from the edge is 
<~ (4 — 39) x 10 23 cm~ 2 (assuming a spherical 
distribution of the absorber and a Fe cosmic abun- 
dance relative to hydrogen of ^ 3 xl0~ 5 ) which 
is consistent with ~ (4 - 7) x 10 23 cm~ 2 mea- 
sured assuming an intrinsic origin for the absorber 
(Table 6a). The predicted Fe Ka line equivalent 
width is in the range of ~ 40-350 eV (Makishima 
1986, Leahy & Creighton 1993) consistent with the 
computed upper limits. 

3.3. ROSAT and ASCA Combined Temporal 
and Spectral Analysis 

For five of the nine quasars: S5 0014+81, NRAO 
140, PKS 0438-436, PKS 0537-286 and PKS 2126- 
158, the best-fit spectral parameters derived from 
ROSAT and ASCA data are consistent with each 
other at a 90% confidence level. Despite the fact 
that the ASCA spectra require systematically more 
absorption than the R OS A T spectra (see Appendix 
B), the 0.1-2 keV flux extrapolated from the ASCA 
spectra (Table 5) is in good (within ~ 10%) agree- 
ment with the measured ROSAT flux (Table 4). 
Therefore, the data from the two instruments were 
combined to benefit of the higher sensitivity of the 
ROSAT PSPC at lower energies and of the ASCA 
instrument higher energy resolution and broad band. 
The normalizations of the 2 instruments were left 
free to vary independently. The unprecedented 
quality (for high-z quasars) of these combined spec- 
tra, covering typically an energy range of ~ 0.4 
40 keV in the quasar frame, provide an excellent 
opportunity to test the data against more complex 
emission models such as complex absorption, ther- 
mal emission and reflection models, as described in 
the following. 

3.3.1. Single Power-Law Fits 

Again, we first fitted the combined spectra with 
a single absorbed power law model, with absorp- 
tion abundances and cross-sections from Morrison 
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& Mc Cammon (1983). The results from these 
fits are given in Table 6a. Spectra, residuals and 
contour plots are shown in Fig. 2. The weighted 
mean of the photon indices for the 5 quasars is 
< r .i-iokcV > — 1-62 ± 0.02, with a dispersion 
co.i-iokcV — 0.11. Note that, when compared 
with ASCA results alone (Table 5), the addition of 
ROSAT data for S5 0014+81 and NRAO 140 gives 
almost unchanged results. However, for the other 
three sources: PKS 0438-436, PKS 0537-286 and 
PKS 2126-158, there is indeed a significant improve- 
ment combining the data. Note that the column 
density for PKS 0537-286 turns out to be marginally 
(> 90%) higher than the Galactic value since sta- 
tistical errors have been reduced, suggesting excess 
absorption in this source too. However, considering 
the ASCA SIS systematic uncertainties (Appendix 
B), conclusions on this issue are unwarranted. 

The inclusion of an Fe K emission line or Fe K 
absorption edge in the model gave results almost 
identical to the one obtained in §3.2 and given in 
Table 5. 

3.3.2. Absorption Fits 

Given sufficient energy resolution, it should in 
principle be possible to constrain the redshift and 
metal abundances of a neutral absorber through di- 
rect spectral fitting. Though such measurements 
sound difficult with present generation X-ray tele- 
scopes, it might be worth trying them with the best 
available data. 

As a first step, therefore, elemental abundances 
were fixed at the cosmic value, given by Morrison & 
Mc Cammon (1983), and spectral fittings with two 
separate absorbers, one at z=0 fixed at the Galactic 
value, and one at the quasar redshift with Ah free 
were repeated. The resulting parameters are given 
in Table 6a. Accordingly, the column densities ob- 
tained for an absorber intrinsic to the quasar are 
larger, ranging between ~ 1.7-55 x 10 21 cm~ 2 . If 
the redshift of the absorber is left free to vary, no 
preferential solution is found for any of the quasars 
indicating that the present data do not allow to dis- 
tinguish between a local absorber (at z=0) or any 
other absorber placed between us and the quasars. 
Fixing the absorber at the quasars redshifts, if the 
abundances of the metals (C, N, O, Ne, Mg, Si, 
S, Ar, Ca, and Fe) were tied together and left free 
to vary (increasing by one the number of degrees 
of freedom) , no relevant constrain was obtained on 
the abundances. This indicates that the model was 
too complex for the present data and no informa- 
tion on the metal abundances of the absorber could 
be obtained. 

The spectra were then fitted again with a single 
absorbed power law model but this time the Mor- 



rison and McCammon cosmic abundances for the 
photo-electric absorption (placed at z=0) were re- 
placed with more recent abundances values obtained 
by Anders & Grevesse (1989) and Feldman (1992). 
The two cases gave systematically lower absorption 
column densities by up to ~ 12% and 10%, respec- 
tively. Since these differences were, in every case, 
within the reported statistical errors, they do not 
substantially modify the present conclusions. This 
systematic difference should, however, be kept in 
mind in particular for future, more sensitive, X-ray 
missions. 

It should be pointed out, however, that the ab- 
sorber could be ionized (Krolik & Kalhnan 1984) 
as is commonly observed in Seyfert 1 galaxies (e.g. 
Otani 1995). If present, a warm absorber should 
imprint two main characteristic features in the low 
energy spectrum: absorption edges at ~ 0.7-1 keV 
(quasar frame) due to OVII and OVIII and extra 
emission below ~ 0.7 keV due to the reduced warm 
absorber opacity at low energies. None of these 
features are observed in the present data. How- 
ever, this is not surprising since, for a quasar at 
z = 2, the energy range for the edge features is 
redshifted to ~ 0.2-0.4 keV which falls outside the 
ASCA energy range and the extra emission is ex- 
pected below ~ 0.2 keV where the PSPC response 
drops down rapidly. Moreover, since the present 
spectra indicate that the ROSAT spectra were typ- 
ically less absorbed than the ASCA spectra (§3.2), 
this may be somehow the signature of a warm ab- 
sorber. In order to test this hypothesis, a warm 
absorber model ("absori" in the XSPEC package) 
placed at the quasars' redshift has, therefore, been 
directly fitted to the data. Free parameters were 
the ionization parameter (£ = L/nR 2 (ergs s _1 cm), 
where n is the warm gas density and R the distance 
of the absorber from the source), and the warm col- 
umn density Ny? (cm -2 ). A temperature of 3 x 10 4 
K (Reynolds & Fabian 1995) and Fe cosmic abun- 
dances were assumed. Results are reported in Table 
6b. Though the fits are generally statistically worse 
than with a single absorbed power law model, the 
data accumulated to date do not rule out this model 
and are consistent with ionization parameters of the 
order of ~ few tens up to few hundreds ergs s -1 cm 
with warm column densities ranging between ~ 2-7 
x 10 22 cm -2 . 

3.3.3. Alternative continuum models 

The possibility that other continuum models could 
be applied to the data was considered in the attempt 
to explain the low energy cut-off with Galactic ab- 
sorption only and/or provide an alternative descrip- 
tion of the higher energy power law spectrum. 

A model consisting of a broken power law with 
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the absorption column density fixed at the Galac- 
tic value was first applied in order to check whether 
a flat power law is able to describe the observed 
low energy cut-off. The spectral curvature may, in- 
deed, be a property of the intrinsic emission of the 
quasars, as a break in the continuum is expected 
if the emission is due to the synchrotron mecha- 
nism where radiative losses are likely to steepen the 
spectra at higher energies. Similar arguments are 
often used to explain the convex shape of the X- 
ray spectra of BL Lac objects (e.g. Barr, Giommi 
& Maccagni 1988). The results of these fits are 
shown in Table 6b. In order to obtain acceptable 
fits, the soft component is required to be very flat 
with typical values of the photon index < 1, while 
the photon index of the hard component is basically 
identical to the values found with a single power 
law model. In the four cases where the break en- 
ergies are constrained, they are found to be in the 
~ 0.7-1.5 keV energy range in the observer frame. 
These fits are statistically acceptable and, in princi- 
ple, very flat spectra could be explained, for exam- 
ple, in the framework of inverse Compton models as- 
suming a particular energy distribution of the elec- 
trons population (i.e. with a convex shape produced 
by the effect of radiative cooling at high energy and 
escape at low energies (Ghiscllini, private commu- 
nication)). This model therefore cannot be directly 
ruled out. However, considering that the broken 
power law model requires an extra free-parameter 
and that, statistically, it is never significantly bet- 
ter than the single absorbed power law model, using 
this model is unjustified with the present data. 

The data were then fitted with a thermal brems- 
strahlung model. With the absorption fixed at the 
Galactic value, this model is also not acceptable be- 
cause there are significant and systematic deviations 
in the residuals at low (E < 1 keV, observer frame) 
energies. With free absorption, we find that the 
spectra can also be described by a very high tem- 
perature (kT ~ 20-90 keV in the source frame) ther- 
mal model and an absorption slightly reduced (but 
still significantly (>99%) higher than the Galac- 
tic column) when compared to the single absorbed 
power law model (Table 6c). The fits are, how- 
ever, in all cases worse than with a single power 
law model. Note, moreover, that at these tempera- 
tures, a bremsstrahlung emission model is virtually 
indistinguishable from a flat power law in the ASCA 
energy range and as already noted by S96, these 
temperatures correspond to ~ 10-15 keV in the ob- 
server frame, which is suspiciously similar to the 
ASCA higher energy limit. Since the temperature 
of bremsstrahlung emission is primarily determined 
by its high energy cut-off, only relatively poor con- 
straints can be set with ASCA on this parameter 



(errors of ~ 25-90%, see Table 6c). As a result, 
the description given with the power law model is 
preferable, and the thermal model will not be con- 
sidered any further. 

Previous Ginga results (Nandra & Pounds 1994) 
have shown that the canonical X-ray slope (r ~ 1.7) 
of Seyfert 1 galaxies can be interpreted as the sum 
of an intrinsic steep power law with r 2 -20kcV ~ 
1.8 — 2.1 plus a reflection component (e.g. Light- 
man & White 1988). Prompted by these results 
and the fact that the present sample of quasars ex- 
hibit fairly flat (r ~ 1.5 — 1.7) spectra up to ~ 
30-40 keV (rest frame) , we searched for evidence of 
high-energy flattening. A reflection component was 
included into the absorbed power law model leaving 
the absorption and photon index free to vary. In its 
simplest form, this model (plrefl in XSPEC) adds 
only one free parameter: the relative normalization, 
R (=A re fl/A p i), between the reflected component 
(A re fl) and the incident power law (A p i) (see Cappi 
et al. 1996 for more details on this model). The 
reflected component was integrated over all view- 
ing angles. Present results give no evidence for a 
reflection component, with upper limits of R rang- 
ing from ~ 0.3 to ~ 0.7, at 90% confidence limit 
for two interesting parameters (Table 6c). The lack 
of Fe K emission lines further supports this conclu- 
sion though it should be reminded that broad lines 
cannot be excluded from the data. 

In order to test the consistency of the present 
results with what is commonly found in Seyfert 1 
galaxies, we forced the intrinsic power law to be 
steep, with a photon index fixed at F = 1.9 and 
fitted the data again. As a result, the model spectra 
required a substantial amount of reflection, R ~ 1.1 
to 3.1, in order to explain the observed flat spectra. 
However, since all the fits became significantly worse 
(Ax 2 ~ 10-20), the presence of such a reflection 
component can be ruled out on statistical basis. 

On the other hand, a high energy cutoff has been 
discovered with GRO OSSE in the intrinsic spec- 
trum of Seyfert galaxies which can be described 
by an exponential law of the form E _r exp(-E/E c ) 
(Jourdain et al. 1992) but where the precise value 
of the e- folding energy, E c , is still uncertain (~ 40- 
300 keV). A hard X-ray cutoff was therefore added 
to the model and fitted to the data but no con- 
strain was obtained for the e-folding energy, with 
lower limits on E c ranging between 15-40 keV in 
the quasar frame. 

4. Comments on Individual Objects 
4.1. S5 0014+81 

S5 0014+81 is the farthest and optically bright- 
est quasar in this sample. It was observed in X- 
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rays by EXOSAT in 1984 (Lawson et al. 1992), 
by the ROSAT PSPC in 1991 (E94, Bechtold et 
al. 1994) and by ASCA in 1993 and 1994. The 
GIS data obtained from the first of the two ASCA 
observations are discussed by Elvis et al. (1994c). 
Their results are in very good agreement with the 
present GIS results. The ROSAT image indicates 
that there are two X-ray sources within ~ 6' of the 
quasar. One source at about 1.8' from the quasar 
was excluded from the ROSAT analysis and was es- 
timated to be negligible in the ASCA energy band. 
The other source, a V — 8.8 KO star at ~ 5' from 
the quasar was neglected for both instruments. The 
photon index and column density derived here (Ta- 
ble 5) are consistent with previous measurements. 
Considering the instrumental uncertainties and dif- 
ferent instrumental band-passes, the EXOSAT and 
-R05.AT PSPC spectra are consistent in both spec- 
tral index and flux to the ASCA results. However, 
it is the first time that the column density is clearly 
constrained at a value higher than the Galactic col- 
umn at a high (>99.99%) significance level (Fig. 2). 
The absolute amount of excess absorption is ~ 13.4 
x 10 20 cur 2 . It should be noted that S5 0014+81 
is well-known for its Lyman Limit absorbers which 
have been object of extensive studies (e.g. Steigman 
1994 and references therein), but that no damped 
Lya system is known along the line of sight to this 
quasar (Lanzetta, Wolfe & Turnshek 1995). 

4.2. PKS 0332-403 

This object was observed in the ROSAT All-Sky 
Survey (RASS) with a flux between 0.1-2.4 keV 
of 1.75 ± 0.6 x 10~ 12 erg cm~ 2 s _1 (Brinkmann, 
Siebert & Boiler 1993). This is in reasonable agree- 
ment with the ASCA flux of ~ 0.96 ± 0.1 x 10~ 12 
erg cm~ 2 s _1 extrapolated between 0.1-2.4 keV and 
corrected for absorption. The ASCA observations 
were studied by S96, who derived spectral parame- 
ters consistent with our results. 

4.3. NRAO 140 

NRAO 140 has a long history of X-ray observa- 
tions. Previous HEAO-1, Einstein, EXOSAT and 
Ginga observations have shown that its X-ray spec- 
trum is well described by a single absorbed power 
law with L - 1.5-1.8 (Marscher 1988, Ohashi et 
al. 1992). This source is known for showing a large 
and variable X-ray absorption column (Ah ~ 3— 
20 x 10 21 cm~ 2 ) which has been interpreted as the 
passage of Galactic dense clouds across the line of 
sight (Bania, Marscher & Barvainis 1991). ROSAT 
and ASCA observations are discussed in Turner et 
al. (1995). Their combined fit yields L ~ 1.73 
± 0.03 and A H ~ 3 ± 2 x 10 21 cur 2 which is 
nearly identical to the present results (Table 6A). 



The X-ray absorption column is significantly larger 
than the Galactic HI column density (Ari — 14.2 
x 10 20 cm~ 2 ; Elvis et al. 1989) derived from mea- 
surements at 21 cm. However, 12 CO emission mea- 
surement toward NRAO 140 implies a molecular hy- 
drogen column density of (Ah 2 — 17 x 10 20 cm~ 2 ) 
which, when added to the HI column density, is in 
excellent agreement with the observed X-ray col- 
umn (Bania et al. 1991; Turner et al. 1995). These 
results confirm that the excess and variable absorp- 
tion in NRAO 140 plausibly originates from Galac- 
tic molecular clouds passing accross the quasar line 
of sight. It should be noted, however, that these 
results are based on the assumption that the CO- 
to-H2 conversion factor is ~ 3 x 10 20 cm~ 2 K _1 
km -1 s, as commonly found for clouds in the Galac- 
tic plane. But this value is rather uncertain in the 
case of high-latitude molecular clouds such as those 
toward NRAO 140. Assuming the much lower con- 
version factor of ~ 0.5 x 10 20 cm~ 2 K _1 km -1 s ob- 
tained by de Vries, Heithausen & Thaddeus (1987) 
or Heithausen et al. (1993), the molecular hydrogen 
column density is reduced by a factor of ~ 6 and 
implies, again, excess absorption. It should also be 
noted that absorption intrinsic to the quasar is not 
ruled out by the current data. 

4.4. PKS 0438-436 

The ROSAT observations of this object were dis- 
cussed in great detail by Wilkes et al. (1992) and 
E94 and its spectral energy distribution is shown in 
Bechtold et al. (1994). As E94 noted, the apparent 
flux decrease by ~ 30% between the two ROSAT 
observations (Table 4) can be explained if one takes 
into account the different wobble mode between the 
observations. Preliminary results of the ASCA ob- 
servation are given in Serlemitsos et al. (1994) and 
are consistent with those presented here. 

4.5. PKS 0537-286 

The ROSAT and ASCA observations of PKS 
0537-286 were previously discussed by Buhler et al. 
(1995) and S96, yielding results consistent with the 
present analysis. As shown in §3.3.1, combining 
the ROSAT and ASCA spectra, there is marginal 
evidence for excess absorption also in this quasar, 
though, within the uncertainties in the low-energy 
response of the SIS. Note, however, that even after 
removing all the SIS data below 1.5 keV, the excess 
absorption remains statistically significant at > 90% 
confidence level. 

4.6. S5 0836+71 

As Brunner et al. (1994) noted, S5 0836+71 un- 
derwent a flux decrease by a factor of ~ 2 between 
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the two ROSAT observations (Table 4), with no 
evidence of spectral variation. The time scale of 
the variation corresponds to ~ 0.2 yr in the quasar 
frame. The present analysis shows that the ASCA 
flux is consistent with the ROSAT flux in the lower 
state. Compared with the ROSAT spectra, the 
ASCA results provide strong evidence for variable 
absorption (AAT H ~ 8 x 10 20 cm~ 2 ) in the direc- 
tion of S5 0836+71 on a time-scale of less than 2.6 
yr. Including in the model the Galactic absorption 
and assuming that the extra absorption is intrin- 
sic to the quasar, the ASCA fit gives an intrinsic 

Nb. ^ l-18±o'.37 x 1022 cm ~ 2 - Thc change in 
absorption corresponds then to a A7Vh ~ 1 x 10 22 
cm~ 2 on a time-scale of ^ 0.8 yr in the quasar frame. 
It should be emphasized that fitting the ASCA spec- 
tra with a broken power law model absorbed by the 
Galactic column gives r so f t < 0.65, Eb re ak 0.96^q'"o 
keV and Thard — l-37+ ; 03 , with a fit formally ac- 
ceptable (x 2 ed /dof ~ 0.87/502). However, as in thc 
other cases where this model was applied (§3.3.3), 
this model is unlikely (but not ruled out) since it 
requires the soft component to be extremely flat. 
The ROSAT and ASCA spectra were also fitted si- 
multaneously with a single warm absorber model 
(§3.3.2). The fits were, however, not acceptable be- 
cause of large residuals in the ASCA data between 
0.5-1.5 keV. Even a complex model could not, there- 
fore, explain the spectra from the two instruments 
contemporary. Variable absorption (neutral or ion- 
ized) is therefore required by the present data. It is 
interesting to note that S5 0836+71 has also been 
observed several times with EGRET in the high en- 
ergy 7-ray band (Nolan et al. 1996). The reported 
spectrum is very steep (r ~ 2.5 ± 0.5) regardless of 
the (variable) 7-ray intensity. 

4.7. PKS 1614+051 

PKS 1614+051 was marginally detected with the 
Einstein IPC (Wilkes et al. 1994) and in the RASS 
(S96). The source is very faint also in thc ASCA 
observation with a 2-10 keV flux of ~ 2 x 10~ 13 
erg cm~ 2 s _1 . 

4.8. PKS 2126-158 

PKS 2126-158 was first detected in X-rays with 
the Einstein IPC (Zamorani et al. 1981). The 
ROSAT spectrum is discussed in E94 and the spec- 
tral energy distribution is given in Bechtold et al. 
(1994). Note that two additional ROSAT observa- 
tions have been analyzed in the present work which 
almost doubled the total number of ROSAT counts. 
The ASCA observations were discussed by Serlcmit- 
sos et al. (1994), who tentatively constrained the 
redshift of the absorber at z<0.4. However, we can- 
not reproduce such contours, neither with nor with- 



out the addition of thc ROSAT data,. No redshift is 
preferred from our analysis. This discrepancy might 
be attributed to the fact that Serlemitsos et al. use 
older response matrices obtained from the prelimi- 
nary calibrations of the GIS and SIS instruments. 

4.9. PKS 2149-306 

RASS and ASCA observations of PKS 2149-307 
were discussed in S96. The comparison of the spec- 
tra obtained from the two instruments indicates 
spectral variability in this source, most likely in- 
terpreted in terms of variable absorption (see Fig. 
2 of S96). The column density found in the present 
analysis is consistent with, but slightly higher than, 
the value obtained by S96. However, it should 
be noted that, unlike the S5 0836+71 case, the 
extra-absorption is not very large in this source 
(AJV H - 4 x 10 20 cm" 2 in S96, AA H - 6 x 10 20 
cm -2 in the present analysis) and could, therefore, 
be affected significantly by the SIS calibration un- 
certainties discussed in §3.2 and Appendix B. 

5. Discussion 

5.1. Excess Absorption 

The most striking result obtained from the present 
analysis is that the X-ray spectra of 6 high-z RLQs 
examined here are absorbed by column densities sig- 
nificantly higher than the Galactic value. Moreover, 
since these measurements correspond to the higher 
signal-to-noise data, they suggest that absorption 
is common (maybe ubiquitous) in high-z RLQs. In 
this section, we discuss the possible origin for such 
absorption. 

5.1.1. Line- of- Sight Absorption 
a) Galactic 

Assuming that the extra absorption originates 
somewhere between the observer and the quasars, 
the first thing to examine is whether all of it could 
be produced in our Galaxy. This evaluation is com- 
monly done through radio measurements at 21 cm 
which have been found to reasonably well trace the 
total amount of gas in the Galaxy (Dickey & Lock- 
man 1990). However, the column densities are de- 
rived from measurements averaged over fairly large 
areas (several tens of arcmin 2 ). Moreover, radio ob- 
servations only detect interstellar atomic hydrogen 
and neither molecular gas (e.g. H 2 , CO, OH, etc.) 
nor dust. So it may happen that, for a given path, 
the total effective absorbing column is significantly 
higher than that indicated by 21 cm measurements. 
If one also considers the proper motion of molecular 
clouds, it could well be possible to explain strong 
and variable absorption by means of special condi- 



9 



tions in our Galaxy alone. A remarkable example is 
the case of NRAO 140 which is located at the edge of 
a molecular cloud (the IC 348 cloud in the Perseus- 
Taurus region, Ungerechts & Thaddeus 1987). As 
discussed above (§4.3), when the molecular content 
along the line of sight is properly taken into account 
through 12 CO measurements to calculate the total 
effective absorption column, it is possible to account 
for the large X-ray absorption column measured for 
this quasar. Under this hypothesis, the absorption 
variability may plausibly be due to the passage of 
absorbing clouds across the line-of-sight. Given this 
possibility, we checked for Galactic CO emission in 
the direction of all quasars of the present sample 
(Table 7). Surveys of local CO emission toward 
extra-galactic sources did not detect CO emission 
in the direction of S5 0836+71 and PKS 2149-306 
(Liszt & Wilson 1993, Liszt 1994). A deep (rms ~ 
0.04 K) observation in the direction of PKS 2126- 
158 puts another limit on any contribution from 
Galactic absorption in this source (Hartmann, pri- 
vate communication). From Table 7, note that S5 
0014+81 is the source with the lowest Galactic lat- 
itude after NRAO 140 which increases the proba- 
bility for a contamination from Galactic molecular 
clouds. Furthermore, S5 0014+81 is located at the 
edge of a large molecular cloud complex (the Polaris 
Flare) and CO emission is detected from positions 
nearby the quasar (e.g. from a point of local maxi- 
mum at 1=120.50 and b=18.63) over an area of more 
than 1 deg 2 (Heithausen et al. 1993). It is difficult 
to precisely quantify the contribution of the cloud 
to the total column density along the line of sight 
to S5 0014+81. However, we believe it is probably 
negligible in this case because from the CO emission 
map shown in Fig. 5 of Heithausen et al. (1993), 
the CO line intensity should be 5; 0.4 K km s _1 in 
the direction of S5 0014+81. A rather conservative 
assumption for the CO-to-H2 conversion factor of ~ 
2-3 x 10 20 molecules cm _2 K _1 km _1 s (Strong et al. 
1988) would yield, therefore, a column density Nn 2 
~ 1.6 x 10 20 atoms cm , which is in any case neg- 
ligible if compared to the atomic hydrogen column 
density of N H ^ 13.9 x 10 20 atoms cm~ 2 derived 
from 21 cm observations (Dickey & Lockman 1990). 

Another indirect tracer of the total column of 
gas along a given path is the emission in the in- 
frared band, in particular at 100 /im. The IRAS 
survey has shown that the Galactic IR emission is 
composed of a diffuse background emission plus sev- 
eral large, filamentary "cirrus" features (Low et al. 
1984) predominantly associated to dust in clouds 
with column densities of a few x 10 20 cm~ 2 , and 
located even at high Galactic latitudes. Using the 
IRAS Faint Source Survey catalog provided on-line 
(Moshir 1989) the IR emission within a radius of 6' 



of each quasar was investigated (Table 7). There 
is no evidence for contamination from cirrus clouds 
in these quasars, except for PKS 2126-158 for which 
the IRAS Faint Source Survey maps indicate a prob- 
able contamination (Wheelock et al 1994). From 
the IRAS maps, the 100 /j,m flux in the direction of 
PKS 2126-158 is estimated to be < 1.9 Jy (the max- 
imum value in that region) which corresponds to a 
brightness < 1.5 MJy sr _1 for a typical IRAS source 
covering 3' x 5'. Adopting a conservative range of 
dSioo/dA H ~ 0.5-2.0 MJy sr^/lO 20 cm" 2 (Reach , 
Heiles, & Koo 1993, Heiles, Reach, & Koo 1988, de 
Vries et al. 1987) as typical dust-to-gas ratios, this 
implies a total (H 2 + HI) hydrogen column density 
of 0.7-3 x 10 20 cm~ 2 . These values are close to the 
21 cm estimate of A H ~ 4.85 x 10 20 cur 2 (Elvis et 
al. 1989), thus ensure that there is no significant 
excess IR emission and therefore that cirrus con- 
tamination is likely to be negligible for the present 
discussion even in this case. 

In summary, on the basis of the above local CO 
properties and IR emission in the direction of these 
quasars, it is unlikely that all the extra absorption 
measured in X-rays is attributable to absorption 
by molecular clouds in our Galaxy, except possibly 
for NRAO 140. However, some contamination from 
molecular clouds and cirrus clouds may be present 
in S5 0014+81 and PKS 2126-158 respectively, but 
at a low level. 

b) Extra-galactic 

It has been known for decades that the space den- 
sity of quasars decreases above z~3 (see Shaver et 
al. 1995 for a review). Though one plausible expla- 
nation for this decrease may be related to the way 
quasars evolve, it has also been argued that dust and 
gas in intervening galaxies could explain or at least 
contribute to the apparent turnover in space density 
(Ostriker & Heisler 1984). Examples of probable in- 
tervening galaxies have already been found in sev- 
eral high-z BL Lac objects (PKS 1413+135: Stockc 
et al. 1992; AO 0235+164: Madejski 1994; PKS 
0537-44 and Wl 0846+561: Narayan & Schneider 
1990; MS 0205.7+3509: Stocke, Wurtz & Perlman 
1995). Interestingly, soft X-ray absorption is often 
present at a level of ~ few x 10 21 cm -2 , which is 
similar to our findings. 

A related possibility is absorption by gas and/or 
dust in damped Lya systems (e.g. Fall & Pei 1993, 
1995) which are plausibly associated to galaxy pro- 
genitors along the line of sight (Wolfe 1995). As ex- 
tensively discussed by Elvis and collaborators in a 
variety of papers (E94 and references therein) which 
included three of the present quasars (PKS 0438- 
436, PKS 2126-158 and S5 0014+81), the numbers 
and column densities of intervening damped Lya 
systems (Lanzetta et al. 1995) could possibly ex- 
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plain the measured X-ray absorption. It should be 
noted, however, that if the hypothesis of the inter- 
vening absorber is correct, whether it is a galaxy or 
a damped Lya system, it predicts that the absorp- 
tion should remain constant over long time-scales. 
This is in apparent contrast with our findings where 
significant absorption variability has been detected 
in S5 0836+71 (§4.6) though this "blazar-type" ob- 
ject may be intrinsically different from the other 
quasars. Moreover, an absorption edge has been 
marginally detected in S5 0014+81 at the energy 
expected if the absorber is at the quasar redshift. 
These results suggest that the absorption is more 
likely intrinsic to the quasars, as discussed in the 
following. 

5.1.2. Intrinsic Absorption 

a) Absorption by dust and/or neutral gas 
Recent observations are reviving the long-standing 
idea (Rieke, Lebofsky & Kinman 1979, Sanders et 
al. 1989) that quasars themselves may be embedded 
in large quantities of dust or thick gas. Webster et 
al. (1995) have recently claimed that a large popu- 
lation of radio-loud quasars are so red that they may 
have been missed by optical searches. This redden- 
ing is interpreted as arising from dust which is most 
likely located within the quasar host galaxy since 
the effect appears to be independent of redshift. 
Other evidence for obscuration have been found in 
several type 2 AGNs, including several narrow line 
radio galaxies of low and high luminosities (e.g. Cen 
A: Bailey et al. 1986; 3C 109: Goodrich & Cohen 
1992, Allen & Fabian 1992; Cygnus A: Ucno et al. 
1994; 3C 265: Dey and Spinrad 1996). In these ob- 
jects, the absorption is thought to arise in a nuclear 
torus, with a wide range of column densities (up to 
./Vh ~ 10 24 ~ 25 cm -2 ). These arguments argue in 
favor of the presence of large quantities of obscur- 
ing material, possibly in the form of dust. Other 
indirect evidence for large quantities of matter in 
high-z objects comes from the recent measurement 
by Ohta et al. (1996) of a strong CO (5-4) emission 
line from a z=4.69 radio-loud quasar. Assuming a 
Galactic CO-to-H 2 ratio, they infer a total mass of 
cold gas ~ 10 11 Mq. 

Another possibility for the site of neutral ab- 
sorbing material is in a cluster of galaxies. Since 
RLQs at z>0.5 are sometimes located in rich clus- 
ters of galaxies (Yee & Ellingson 1993), cluster cool- 
ing flows may provide the necessary cold material 
to produce the observed X-ray absorption (White 
et al 1991, E94). The radio galaxy Cygnus-A is 
a remarkable example of such scenario (Reynolds 
& Fabian 1996). However, again, if the variable 
absorption observed in S5 0836+71, and possibly 
PKS 2149-306 is real, this could not be explained 



either by an absorption torus nor by a cooling flow 
model. Thus an alternative explanation must be 
found, at least for these cases. One possible expla- 
nation could be the passage of neutral material (e.g. 
dense absorbing clouds) across the radiation field 
(possibly anisotropic). On one hand, cold clouds or 
"blobs" are expected to survive the intense radia- 
tion field close to the central regions of AGNs (Guil- 
bert & Rees 1988, Cclotti, Fabian & Rees 1992) and 
may well imprint reflection features (iron line and 
hard tail) commonly seen in the X-ray spectra of 
AGNs (Nandra & George 1994). On the other hand, 
anisotropic radiation, possibly with beaming of the 
radiation, may be expected in these objects on the 
basis of the extremely high luminosity (L2-iokcV ~ 
10 47-48 erg s _1 , in the quasar frame) observed and 
on the basis of general arguments on the unification 
of radio-loud quasars (Urry & Padovani 1995). As 
far as S5 0836+71 is concerned, supcrluminal ra- 
dio components (Krichbaum et al. 1990), a rapid 
7-ray flare observed by EGRET (von Linde et al. 
1993), and theoretical arguments (Dondi & Ghis- 
ellini 1995) support the hypothesis that, at least in 
this object, the radiation source is relativistically 
beamed. Optically thin clouds with column densi- 
ties of JVh ~ io 21 ~ 23 cm~ 2 covering the continuum 
radiation or, alternatively, optically thick (Ah > 
10 25 cm~ 2 ) clouds covering quite a large fraction 
of it could plausibly reproduce the observed X-ray 
absorption. The variability could arise from the pas- 
sage of the clouds across the radiation field line of 
sight. 

b) Absorption by ionized gas 

The absorption may also be produced by par- 
tially ionized "warm" material. Warm absorbers 
have been previously invoked to explain the complex 
X-ray absorption observed in several active galaxies 
(Nandra & Pounds 1994; Cappi et al. 1996) and 
quasars (MR 2251-178: Halpern 1984, Otani 1996; 
3C 351: Fiorc et al. 1993, Mathur et al. 1994a; 3C 
212: Elvis et al. 1994d, Mathur 1994b). The state 
of a warm absorber depends only, though critically, 
on the ionization parameter, £ = L/nR 2 (Kallman & 
McCray 1982). Since warm absorbers have been un- 
ambiguously detected in low luminosity AGNs, they 
may survive also at higher luminosities since, for a 
given density, the gas may just be at a larger dis- 
tance from the central source. As shown in §3.3.2, 
the present data are consistent with a warm ab- 
sorber model as a possible explanation of the low 
energy absorption. However, the quality of the data 
coupled with the high redshifts do not allow us to 
distinguish between a cold and ionized absorber. It 
is interesting to note that if the ASCA X-ray spec- 
trum of MR 2251-178 (z ~ 0.068; L 2 _i 0k cV ~ 10 45 
erg s^ 1 ), which clearly indicates the presence of a 
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warm aborber (£ ~ 28, iV warm ~ 2.6 x 10 21 cm -2 ), 
is redshifted to z=2 (i.e. the data are cut at E > 
1.2 keV), then it could be equally well fitted with 
a neutral absorption column density of <~ 6-12 x 
10 20 cm -2 (Otani, private communication), similar 
to what we find in our sample. This implies that 
the data accumulated to date for high-z quasars can 
hardly be used to distinguish between a cold or ion- 
ized absorber. A warm absorber cannot, therefore, 
be ruled out. 

If the warm absorber hypothesis is correct, vari- 
able absorption could be interpreted as variation in 
the ionization state of the absorber (e.g., MCG-6- 
30-15: Otani et al. 1996). An ionized absorber also 
predicts an absorption edge at E > 7.1 keV which 
seems to be at variance with our findings for S5 
0014+81 where the fitted energy of the absorption 
edge is centered at ~ 7.1 keV. The quality of the 
data, however, does not exclude mildly ionized ab- 
sorption either. 

5.2. Statistical Properties 

Correlations among the following parameters: 
photon index, column densities, 2-10 keV luminos- 
ity, radio-loudness, and source redshift have been 
searched for in the whole ASCA sample. The anal- 
ysis shows no significant correlation which is unsur- 
prising given the small size of the sample. Though 
not adequate for statistical correlations between 
the X-ray emission parameters, the sample is large 
enough to be representative of the entire class of "X- 
ray selected" high-z RLQs. But because the sources 
were selected among the brightest and farthest ones 
detected by previous satellites, the sample is clearly 
biased toward the most luminous objects in the uni- 
verse. 

Earlier X-ray spectra of RLQs, mostly at low red- 
shift, showed mean photon indices between ~ 1.4 
1.8 both in the soft - 0.1-4 keV (Einstein IPC: 
Wilkes & Elvis 1987; ROSAT PSPC: Brunncr et al. 
1994) and hard - 2-10 keV (EXOSAT: Comastri 
et al. 1992; Ginga: Lawson & Turner 1996) energy 
bands. In some works, significant dispersion around 
the mean has been found (Brunner et al. 1994, Law- 
son & Turner 1996). If real, this indicates that ei- 
ther RLQs have an intrinsic dispersion of slopes or 
the intrinsic X-ray spectrum is more complex than a 
single absorbed power law (e.g. soft-excess or com- 
plex absorption). Previous X-ray spectra of high- 
z RLQs taken with ROSAT (E94, Bechtold et al. 
1994b) and ASCA (Serlemitsos et al. 1994, S96) are 
consistent with an average photon index between ~ 
1.5-1.8 up to 30 keV in the source frame. 

The mean photon index derived for the whole 
sample over the full ASCA energy band is < Po.s-iokcV > 
~ 1.61 ± 0.04 with a dispersion cro.5-iokcV — 0-10 =t 



0.03. This average value is in agreement with most 
of the findings described above for low-z and high- 
z RLQs derived in both the soft and hard band, 
and indicates that there is no X-ray spectral evo- 
lution with redshift or luminosity over a broad ~ 

1- 40 keV energy band (source frame). This is fur- 
ther illustrated in Fig. 3 which shows the X-ray 
spectral slopes observed with ASCA between ~ 0.5- 
10 keV as a function of redshift and 2-10 keV lu- 
minosity (quasar frame) together with the photon 
indices obtained with Ginga for 18 low-z RLQs as 
reported by Lawson & Turner (1996). It is interest- 
ing to note that the dispersion is significantly differ- 
ent from zero also at high-z, though it is somewhat 
smaller than the Ginga value a ~ 0.19lg o4 (Law- 
son & Turner 1996). However, the mean photon 
index obtained with ASCA in the 2-10 keV energy 
band, where the present measurements are not af- 
fected by the absorption, is < r2-iokcV > — 1-53 
± 0.05, with a dispersion <72-iokcV 0.12. The 
intrinsic dispersion is, in this case, consistent with 
zero within la. This result implies that the above 
non-zero dispersion is probably a consequence of the 
heterogeneous absorption measured in the data and 
that the intrinsic distribution of photon indices is 
well characterized by a single slope. 

Excluding S5 0836+71 from the sample because 
of its blazar-like properties, the results are as fol- 
low: < r .5_iokeV > — 1-64 ± 0.06 with <t . 5-iokoV 
~ 0.071^^ and < F 2 -iokcV > ^ 1.57 ± 0.09 with 

02- iokcV ~ 0.14. These are consistent with the 
above results. 

5.3. Implications for the XRB and GRB 

It is by now widely accepted that AGNs could 
supply the bulk of the XRB emission above a few 
keV. Detailed modeling, based on the unified scheme 
and assuming the existence of a large population of 
absorbed AGNs, provide a good fit to the observed 
XRB spectrum up to ~ 100 keV (Madau, Ghisellini 
& Fabian 1994; Comastri et al. 1995). If the ab- 
sorption discovered in some of the objects in our 
sample turns out to be a common property of high- 
z, high-luminosity RQQs this would give a strong 
support to the above models. 

It should be noted that RLQs constitute only a 
small fraction of the quasar population so that their 
contribution to the XRB is of the order of a few 
percent. This point is also enforced by our find- 
ings that the observed X-ray spectra of RLQs have 
r ~ 1.5 — 1.6 up to ~ 40 keV in the source frame, 
steeper than the XRB as illustrated in Fig. 3. But 
the derived mean slope for the present sample has 
interesting implications for the 7-ray background. 
In fact, it has been recently pointed out (Comastri, 
Di Girolamo & Setti 1996) that flat spectrum radio 



12 



quasar (FSRQ) may provide a significant fraction 
(~ 70-80 %) of the 7-ray background above several 
tens of MeV. This result has been achieved assum- 
ing that the X-ray spectra of FSRQ in the X-ray 
band (from ~ 1 keV up to several hundreds of keV) 
are characterized by a slope with T = 1.5. Our 
findings provide, therefore, further support for this 
model suggesting that RLQs are likely to provide 
an important contribution to the hard X-ray - 7- 
ray backgrounds. 

6. Conclusions 

The present analysis of ASCA and ROSAT ob- 
servations of high-z RLQs confirms that excess ab- 
sorption is common in these objects. It is indeed 
detected significantly in the six brightest quasars 
of our sample. Unfortunately, even combining the 
two instruments, it is not possible to constrain the 
redshift of the absorber through direct spectral fit- 
ting. However, new results for S5 0836+71, and 
possibly S5 0014+81, favor the hypothesis that, at 
least in these two cases, the absorption is intrin- 
sic to the quasars. Indeed, a comparison of ASCA 
and ROSAT spectra indicates the absorption has 
varied (AiV H - 8 x 10 20 cur 2 ) in S5 0836+71 on 
a time-scale of ~ 2.4 yrs (observer frame) and an 
Fe K absorption edge is marginally detected in the 
ASCA observations of S5 0014+81, at the energy 
expected if the absorber is at the quasar's redshift. 
If we assume that the absorption is intrinsic for all 
objects in the sample, the column densities implied 
range between ~ (1.7-55) x 10 21 cm~ 2 . 

Since the data do not allow to constrain unam- 
biguously the redshift of the absorption, several pos- 
sible origins have been considered. A careful search 
for Galactic CO and/or IR (100 /zm) emission in 
the direction of the quasars indicates that the con- 
tribution to the total X-ray absorption column from 
molecular gas and dust present in our Galaxy is neg- 
ligible for the majority of the quasars. It is, however, 
likely to be significant in the case of NRAO 140. 
Extra-galactic absorption from intervening galaxies 
and/or damped Lya systems may be relevant for 
some of the objects, but is not consistent with the 
new findings for S5 0836+71 and S5 0014+81. In- 
trinsic absorption from an absorption torus or intra- 
cluster gas is also unlikely in the light of the ab- 
sorption variability. It is argued, therefore, that 
the preferable explanation is intrinsic absorption by 
cold gas, possibly in the form of clouds, near the 
central source. A warm absorber would be a plau- 
sible candidate as well. 

The accumulated counts allowed us to put, for 
most objects, stringent limits on the presence of 
an Fe K emission line, with upper limits on the 
EW ranging between ~ 40-400 eV in the quasar 



frame. There is, however, more uncertainty about 
the strength of the line if the width is left free to 
vary. Alternative continuum models have also been 
fitted to the data. The upper limits on the intensity 
of Compton reflection range between R ~ 0.1-0.7. 
A broken power law model, tested as an alterna- 
tive explanation for the low energy cut-off, would 
require unreasonably flat (r < 1) slopes at low en- 
ergies, though it cannot be ruled out by the data. It 
is also argued that a thermal bremsstrahlung origin 
for the observed continuum is unlikely. 

The average photon index (< T > ~ 1.5-1.6) 
is consistent with previous measurements at high 
and low-z, therefore suggesting that RLQs do not 
show spectral changes over about two decades in 
energy (<~ 0.4-40 keV) with either redshift or lumi- 
nosity. This average slope, which is steeper than 
the spectrum of the XRB but flatter than the aver- 
age slopes of RQQs and low-luminosity AGNs, also 
suggests that RLQs do not contribute significantly 
to the XRB but, instead, are likely to provide a sig- 
nificant fraction of the 7-ray background flux. 

APPENDICES 

A. GIS Backgrounds 

In the standard way of analyzing ASCA GIS 
spectra, there are two possible choices for the back- 
ground subtraction: either from a source- free part of 
the observation field of view (hereafter "local back- 
ground"), or from the available blank sky observa- 
tions with the same region filter as used for the 
source spectrum (see "The ABC Guide to ASCA 
Data Reduction"). Both methods should, in princi- 
ple, give similar results. However, slight but signif- 
icant differences in the source spectral slopes were 
found in the present analysis and forced us to choose 
a non-standard background region. 

The problem is that background regions extracted 
from the source field of view almost inevitably fall 
at a substantial off-axis distance because sources are 
usually pointed nearly on-axis. But as the distance 
from the X-ray telescope optical axis increases, re- 
duction of the contribution from the cosmic X-ray 
background (due to the vignetting) and softening 
of the non X-ray background (Kubo ct al. 1994) 
change in a complex and position dependent man- 
ner the spectrum of the background, which cannot 
be corrected in the current analysis procedure. On 
the other hand, standard background regions ex- 
tracted from the blank sky files from the same area 
as that used for the source almost always include 
the regions contaminated by the Seyfert 2 galaxy 
NGC 6552 (Ebisawa 1994). Because of the very 
strong Fe K emission line of this source at ~ 6.4 
keV (Fukazawa et al. 1994), the source spectra 
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obtained with blank sky backgrounds tend to be 
steeper than those obtained with local background 
subtractions. As a result, in the analysis of the 
weakest sources (typically those with intensity < 
0.1 cts/s), different backgrounds often yielded sig- 
nificantly different source slopes with deviations up 
to Ar ~ 0.2. Therefore, we adopted a non-standard 
background extracted from the GIS standard blank 
sky files from a region non contaminated by NGC 
6552 and at the same or, if not possible, similar off- 
axis distance of the source region. This choice of 
background has the advantage to avoid the above 
problems of the off-axis position with the local back- 
ground and of the contamination by NGC 6552 with 
the blank sky background. However, it has the dis- 
advantage that the screening criteria of the blank 
sky background never match exactly those of the 
source observation (similarly to any other choice of 
blank sky background). However, we estimate this 
is a minor problem when standard screening crite- 
ria are adopted because the GIS background repro- 
ducibility (in time) is achieved with systematic er- 
rors less than a few percent (Ikebe et al. 1995). 

B. On the ASCA SIS response uncertainties 
at low energies 

Although the results obtained with ASCA and 
ROSAT always agree within their statistical errors 
(except for S5 0836+71), the present analysis indi- 
cates that the ASCA spectra require systematically 
more absorption than the ROSAT spectra (Fig. 1). 
This suggests that the SIS responses may suffer 
of a systematic excess absorption, though ROSAT 
PSPC calibration uncertainties (e.g. non-linear gain 
variations) may play a role as well. This confirms 
previous, more detailed, studies on the SIS calibra- 
tion uncertainties at low energies (e.g. Hayashida et 
al. 1995). Indeed, it is widely recognized that there 
are local features around 0.5-0.6 keV which show 
up, preferentially in strong sources, in the form of 
absorption edges around 0.5 keV and/or emission 
lines at ~ 0.6 keV (e.g. Guainazzi & Piro 1994, 
Cappi et al. 1996, Dotani et al. 1996). The origin 
of the problem is not clear, though probably related 
to the presence of the Oxygen K-edge (produced in 
the CCD dead layer) at E ~ 0.54 keV in the SIS re- 
sponse. Because of this edge, a small offset and/or 
variations of a few eV in the exact response energy 
scale could easily be responsible for such features. 
Though no clear feature is evident in the low energy 
residuals of the quasars spectra (Fig. 2), it is pos- 
sible that it could have affected somehow our mea- 
surements. Other uncertainties concern the abso- 
lute SIS efficiency at low energies which is currently 
being quantified in terms of a systematic error of 
- 3 x 10 20 cm" 2 for SIS0 and - 2 x 10 20 cm" 2 



for SIS1 in the absorption column (Hayashida et al. 
1995, Dotani et al. 1996). This problem is prob- 
ably related to the correct calibration of the SIS 
quantum efficiency, i.e. the thickness of the CCD 
dead layer. It is not clear, however, how much the 
above local features at E ~ 0.5 keV are responsible 
for and/or related to such excess absorption but, in 
principle, both effects (local features and excess ab- 
sorption) should be already taken into account in 
the systematic error of ~ 3 x 10 20 cm~ 2 . 

In light of the above considerations, we have per- 
formed a series of checks in order to establish the 
reliability of the measured excess absorption for all 
ASCA spectra. 

a) Test 1 - Artificial gain shifts: 

Using the "gain" command in XSPEC, we shifted 
the energy scales of the SIS0 and SIS1 responses be- 
tween ± 20 eV, i.e. within the typical systematic 2a 
uncertainties in the energy scale (Otani & Dotani 
1994). The observed changes in the column densi- 
ties were, in every case, lower than 2 x 10 20 cm~ 2 , 
which is smaller than the typical statistical errors 
found in the present work. 

b) Test 2 - Cut of the data below 0.65 keV: 
When, more drastically, all data below 0.65 keV 

were ignored, the fits yielded even larger best-fit ab- 
sorption columns for the quasars NRAO 140, PKS 
0438-436, PKS 0537-286 and PKS 2126-158. For S5 
0014+81 and S5 0836+71, the columns were lowered 
by ~ 1 x 10 20 cm" 2 (~ 5%) and - 2.5 x 10 20 cm~ 2 
(~ 20%), respectively. In every case, the absorp- 
tion columns remained significantly higher than the 
Galactic value at more than 99% confidence level. 

From a) and b), we conclude that the measured 
SIS excess absorption cannot be attributed to the 
local detector features commonly reported around 
0.5-0.6 keV. The problem is more likely related to 
the SIS efficiency as a whole at low energies. 

c) Test 3 - Comparison with the calibration source 
3C 273: 

The SIS response matrices were calibrated on the 
ground and then cross-calibrated with the GIS in- 
flight using the observations of the quasar 3C 273. 
In particular, the response matrices of SIS0 chip 
number 1 and SIS1 chip number 3 (the chips used 
in the present analysis) were calibrated using an ob- 
servation performed in 1993, December 20 th . The 
calibrations were performed taking the slopes from 
the GIS spectra and assuming Galactic absorption. 
3C 273 was observed 8 more times by ASCA. Re- 
sults obtained from the overall observations (Cappi 
& Matsuoka 1996) do not modify the conclusions 
presented below. For the purpose of investigating 
the efficiency of the SIS, only the results from the 
December 20 th 1993 calibration observation are re- 
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ported below. The observation was analyzed using 
the same criteria (same screening criteria, extrac- 
tion regions, etc.) used for the quasars to allow a 
direct comparison. About 82 000 counts/SIS were 
collected in the 0.4-10 keV energy range, for a to- 
tal exposure time of ~ 15 500 s. The spectra ex- 
tracted for SIS0 and SIS1 were then fitted with a 
single absorbed power law, first separately and then 
simultaneously. Confidence contours obtained from 
each fit are shown in Fig. 4. The best-fit column 
density is larger than the Galactic column in both 
SIS0 and SIS1, yielding an excess column density 
of ~ 2 x 10 20 cm~ 2 in the fit with both detectors. 
This value gives a direct estimate of the systematic 
error we must take into account when interpreting 
the results from the present analysis. This estimate 
is in good agreement with the measurements ob- 
tained by Hayashida et al. (1995). Also, it should 
be noted that independent analysis of the SIS data 
extracted from the central regions of Coma cluster 
also show excess absorption of ~ 3 x 10 20 cm~ 2 
(Hashimotodani, private communication). 

Therefore, we conclude that the total systematic 
error of the absorption column density should be 
conservatively smaller than ~ 3-4 x 10 20 cm~ 2 . In 
the analysis presented in §3 and when discussing the 
results, only deviations larger than this value have 
been considered. 
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TABLE 1: The Radio-Loud Quasar Sample 



Ubject 


R.A. a 




~a 
Z 


Alb 

^Hgal 


m v 


.pa 

^GHz 


R L 




(J2000) 


(J 2000) 




(10 cm ) 




(mJy) 




S5 0014+81 


00 


17 


08.5 


+81 35 08.1 


3.38 


13.9 


16.5 


551 


2.78 


PKS 0332-403 


03 


34 


13.6 


-40 08 25.4 


1.44 


1.43 


18.5 


2600 


4.24 


NRAO 140 


03 


36 


30.1 


+32 18 29.3 


1.26 


14.2* 


17.5 


2500 


3.87 


PKS 0438-436 


04 


40 


17.2 


-43 33 08.6 


2.85 


1.47 


18.8 


7580 


4.84 


PKS 0537-286 


05 


39 


54.3 


-28 39 56.2 


3.10 


1.95* 


20.0 


990 


4.39 


S5 0836+71 


08 


41 


24.4 


+70 53 42.2 


2.17 


2.78 


16.5 


2573 


3.45 


PKS 1614+051 


16 


16 


37.5 


+04 59 33.2 


3.21 


4.90 


19.5 


850 


4.15 


PKS 2126-158 


21 


29 


12.2 


-15 38 40.8 


3.27 


4.85* 


17.0 


1240 


3.33 


PKS 2149-306 


21 


51 


55.5 


-30 27 53.7 


2.34 


1.91 


17.9 


1150 


3.66 



a Coordinates, redshift, V magnitude and radio flux at 5 GHz, from Veron-Cetty & Veron (1993). 

6 Galactic absorption from Dickey & Lockman (1990). The values marked with (*) are from Elvis et al. (1989). 

c Radio-Loudness defined as Rl = log(f5GHz/fv)- We used my(0)=3360 Jy (Wamsteker W. 1981) 



TABLE 2 : ROSAT Observation Log 



Object 


Date 


Wobble 


Matrix 


Exposure (s) 


NC a 


Ref." 


S5 0014+81 


03/15/91 


on 


DRM06 


5951 


394 


1,2 


NRAO 140 


08/08/92 


on 


DRM36 


4039 


992 


3 


PKS 0438-436 


02/19/91 


off 


DRM06 


10725 


645 


1,4 


PKS 0438-436 


09/19/92 


on 


DRM36 


10506 


547 


1 


PKS 0537-286 


09/28/92 


on 


DRM36 


9487 


568 


5 


S5 0836+71 


03/23/92 


on 


DRM36 


6993 


5400 


6 


S5 0836+71 


11/02/92 


on 


DRM36 


5026 


2008 


6 


PKS 2126-158 


05/08/91 


on 


DRM06 


3424 


613 


1,2 


PKS 2126-158 


11/12/92 


on 


DRM36 


3968 


729 


1,2 


PKS 2126-158 


04/27/93 


on 


DRM36 


4160 


741 


This work 


PKS 2126-158 


05/17/93 


on 


DRM36 


1610 


321 


This work 



a NC: Net counts. 

b References: (1) E94, (2) Bechtold et al. (1994), (3) Turner et al. (1995), (4) Wilkes et al. (1992), (5) Biihler et al. 
(1995), (6) Brunner et al. (1994) 
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TABLE 3: ASCA Observations Log 



Object 


Date 


Exposure™ (s) 


Count Rate™ 


(xl0~ 2 s- 1 ) 


CCD mode 


Ref." 






CIS 


SIS 


CIS 


SIS 






S5 0014+814 (AOI) 


10/29/93 


39004 


21750 


7.33 


3.88 


4 CCD 


(1) + This work for SIS 


S5 0014+814 (AOII) 


10/07/94 


27494 


20703 


5.14 


5.98 


1 CCD 


This work 


PKS 0332-403 


08/12/94 


16646 


14153 


1.98 


2.65 


1 CCD 


(2) 


NRAO 140 


02/02/94 


33181 


32178 


16.3 


20.7 


1 CCD 


(3) 


PKS 0438-436 


07/13/93 


34499 


26412 


2.77 


3.20 


4 CCD 


(4) 


PKS 0537-286 


03/12/94 


36236 


29052 


2.51 


4.21 


1+2 CCD 


(2) 


S5 0836+71 


03/17/95 


16532 


10540 


24.0 


30.3 


1 CCD 


This work 


PKS 1614+051 


08/02/94 


39695 


32577 


0.88 


0.55 


1 CCD 


(2) 


PKS 2126-158 


05/16/93 


15914 


14291 


10.1 


12.8 


1 CCD 


(4) 


PKS 2149-306 


10/26/94 


19168 


16504 


19.0 


25.5 


1 CCD 


(2) 



a Reported values for the CIS and SIS are averaged over the detectors (GIS2 with GIS3 and SIS0 with SIS1). 
b References: (1) Elvis et al. (1994c), (2) S96, (3) Turner et al. (1995), (4) Scrlemitsos et al. (1994) 

TABLE 4 : ROSAT X-Ray Spectral Fits 
PSPC - Single Power Law (JVh free) 



Object 


Energy Range" 
(keV) 


N H 
(10 20 cm- 2 ) 


r 




Xred/dof 


Fsx 


log Lg X 


S5 0014+81 


~0.4-10 


20.0t? 6 ;!J 


1 89 +2 - 4U 


5.30 


1.37/17 


0.73 


47.8 


NRAO 140 


~0.2-5 


9 r, = + 21.9 


1 - oo -0.81 


20.7 


1.21/15 


2.88 


46.9 


PKS 0438-436 (Feb 91) 

(Sep 92) 

(tot.) 


~0.4-9 

11 
11 


o -i q+20.4 
0.10_ 3 69 

5.07^:2 
6-3i±t:S 


1 71+ 1 - 16 

L - 1 1 -0.48 
+0.49 
1 - d '-0.44 
i fi o+0.36 
1 -°°-0.32 


3.53 
2.51 


0.45/19 
1.28/25 
0.93/46 


0.71 
0.57 


47.3 
47.0 


PKS 0537-286 


~0.4-10 


q i 1+1.84 


i ^o+ - 45 

l.OU_ 44 


2.61 


1.84/7 


0.65 


47.1 


S5 0836+71 (Mar 92) 

(Nov 92) 

(tot.) 


~0.3-7 

11 
11 


o on+0.48 
°-^ a -().46 
o 4 =+0.79 
°-^°-0.80 

O qq + 0.42 


1 ci+0.13 
1 - Ji -0.12 
i rr+0.22 
l.OO_ 22 
i tr 9 +0.12 


35.0 
17.9 


1.27/18 
1.55/18 
1.34/38 


8.59 
4.41 


47.8 
47.5 


PKS 2126-158 (May 91) 

(Nov 92) 

(Apr 93) 

(May 93) 

(tot.) 


~0.4-10 

11 
11 
11 
11 


Q qo + 13.7 

°- J8 -3.82 

Q 01 + 22 S 
y - ul -4.25 

7 9 q+8.27 
' -^ J -2.89 
11 r + 37.6 
11 -°-6.61 

8 42+ 4 ' 88 

°-^-2.30 


1 44+ - 81 
-■-•^-o.so 

1 ^O 4 " 1 ' 22 
r-ou_ 46 

1 =1+0.57 
1 - J1 -0.42 
1 Q1+2.07 
1 - 31 -0.74 
1 =q+0.34 
l.OO_ 25 


10.4 
10.9 
10.0 
12.7 


0.95/18 
0.93/21 
0.85/15 
1.05/15 
0.90/75 


2.08 
2.13 
2.08 
2.23 


47.7 
47.8 
47.7 
48.2 



Note: Intervals are at 90% confidence, for two interesting parameters. 
a Approximate energy range, in the quasar rest frame. 

b Unabsorbed flux at 1 keV (observed frame) in units of 10~ 4 photons cm -2 s -1 keV 
c Absorbed flux between 0.1-2 keV (observed frame) in units of 10 -12 ergs cm~ 2 s -1 
d Intrinsic luminosity between 0.1-2 keV (quasar frame) in units of ergs s _1 . 
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TABLE 5: ASCA X-Ray Spectral Fits 
GIS+SIS - Single Power Law (TVh free; A^ H = AWO 



Object 


Energy Range" 
(keV) 


AT H 
(10 20 cnr 2 ) 


r 


A fc 


Xred/dof 




*SX' 


log L|, x 


EW (FoK) 1 
(eV) 


T cdgc {FeK) 3 


S5 0014+81 (Oct 93) 

(Oct 94) 

(tot.) 


~2-40 

11 
11 


28.4+5 ^ 

26.9 + g : ? 

28.llt.2 
13.9 Fixed 


1 -'°-0.08 

1 66+ 010 

-, 79 +0.06 
± -'^-0.07 

l-54±o'.ra 


6.46 
6.60 


1.03/314 
0.92/237 
1.02/553 
1.13/554 


2.43 
2.49 


0.82 
0.73 


47.8 
47.7 


<75 
<131 
<70 


<0.23 
21+ 018 
15+ 012 

U - 1O -0.11 


PKS 0332-403 


~l-20 


< 22.9 
1.43 Fixed 


-, 7 n+0.28 
1 - 1 J -0.25 

1.60±g;iJ 


2.48 
2.05 


0.77/103 
0.78/104 


0.94 


0.47 


46.3 


<415 


<0.66 


NRAO 140 


~l-20 


31.8+1° 
14.2 Fixed 


i 7 n+0.05 

1.46±2;{g 


21.8 
16.3 


0.96/864 
1.23/865 


9.09 


2.56 


47.1 


<36 


<0.10 


PKS 0438-436 


~2-38 


1.47 Fixed 


1.63+° 

— 0. 14 


3.89 
2.92 


0.93/236 
1.00/237 


1.75 


0.63 


47.3 


<240 h 


<0.29 h 


PKS 0537-286 


~2-40 


1.95 Fixed 


1 46+ ' 14 

1 Qfi+ U - Utl 

1.00_Q Qg 


3.09 
2.71 


0.74/197 
0.76/198 


1.85 


0.63 


47.4 


<139 


<0A5 h 


S5 0836+71 


~l-30 


ii-4i 3 2 :2 

2.78 Fixed 


1 4^+0-05 

i no+0.03 
J -°^— 0.03 


23.3 
19.3 


0.87/503 
0.97/504 


14.0 


4.25 


47.8 


<110 


<0Al h 


PKS 1614+051* 


~2-40 


4.90 Fixed 


1.6 Fixed 


0.43 




0.21* 


0.10* 


46.4* 






PKS 2126-158 


~2-40 


13.6t|;| 
4.85 Fixed 


1 RS+ ' 09 
1-08-0. 09 
1 rq+0.04 
1 -°°-0.04 


13.1 
10.8 


1.01/331 
1.07/332 


5.51 


2.20 


47.9 


<107 


<0.08 


PKS 2149-306 


~l-30 


o 0+I.8 
o.o_2 2 

1.91 Fixed 


1 54+0-05 

i -°*-0.05 
1 -^ z -0.02 


19.2 
16.3 


1.07/633 
1.13/634 


9.94 


3.81 


47.8 


<85 h 


<0.17 h 



Note: Intervals are at 90% confidence, for two interesting parameters when Nu is free and for one interesting parameter 
when Nn is fixed at the Galactic value. 

a Approximate energy range, in the quasar rest frame. 6 SIS unabsorbed flux at 1 keV (observed frame) in units of 
10~ 4 photons cm~ 2 s _1 keV -1 . c Absorbed flux between 2-10 keV (observed frame) in units of 10 -12 ergs cm~ 2 s _1 , 
calculated with the SIS normalization only. d Absorbed flux between 0.1-2 keV (observed frame) in units of 10~ 12 ergs 
cm~ s~ , extrapolated from the best-fit spectra (i.e. with the SIS normalization). e Intrinsic luminosity between 2-10 
keV (quasar frame) in units of ergs s~\ calculated with the SIS only. 1 90% limit on the equivalent width of the Fe 
K line measured in the observer frame, with a — eV and E = 6.4 keV (quasar frame), for one interesting parameter 
(Ax 2 = 2.71). These values were about a factor of two worse if the line widths are as large as 0.5 keV. 9 90% limit on 
the absorption edge depth with energy fixed at 7.1 keV (quasar frame). h Value affected by the calibration uncertainties 
of the SIS between ~ 1.8-2.4 keV. An estimated conservative systematic error of 10 eV (observer frame) for the lines 
and At = 0.05 for the edge depth have been added to these values. 1 Quasar was detected but with too few counts to 
perform a spectral analysis. Fluxes and luminosity are calculated assuming a photon index = 1.6 (i.e. the average value 
from the sample) and Galactic absorption. 
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TABLE 6: ASCA and ROSAT Combined Spectral Fits 
A. Single Power Law (Ah free; Ah at z + A^Hgai ) 



Object 


Energy Range™ 
(keV) 


A H 
(10 20 cm- 2 ) 


r 


Xred/dof 


S5 0014+81 


-0.4-40 


27.4+ 4 ^ 
ooi_ 170 


1 71 +U.Ub 
x - ' 1 -0.07 


1.03/572 
1.03/572 


NRAO 140 


-0.2-20 


3L1+ 2 ? 


1 70+ - 04 

1 - uu -0 .04 


0.97/881 
98/881 


PKS 0438-436 


-0.4-38 


r 01+2.66 
O.81_i. 40 

71 Q+ 44 - 1 


i-4»_ .09 

1 46 +009 


0.97/281 
0.97/281 


PKS 0537-286 


-0.4-40 


9 09+O.97 
1D -'-14.6 


1 oo+0. 08 

±.00_Q Qg 

1 o 7 +0.08 
L - ,:>l -0.08 


0.77/204 
0.77/204 


PKS 2126-158 


-0.4-40 


10.4+|-| 
116+S 


1 63+ 007 

1,UJ -0.08 
1 ro+0.07 

i.oo_ 06 


0.99/411 
1.01/411 



Note: Intervals are at 90% confidence, for two interesting parameters. 
a Approximate energy range, in the quasar rest frame. 

B. Complex models for the absorption 



Warm Absorber (at z) + -/Vkgai Broken Power Law (Ah = AH ga i) 



Object 


A w 
(cm" 2 ) 


r 


t 

(ergs cm 2 s 1 ) 


X^d/dof 


r SO ft 


Ebrcak 

(keV) 


Thard 


Xred/dof 


S5 0014+81 


74« 


1 79+u.ou 


< 191 


1.03/571 


-0.05+^ 


1 04 +u ' is 

1 - u ^-0.12 


1 60 +u(J4 

l.OU_ 04 


1.03/571 


NRAO 140 


305±lf 4 


1 7 9+0.05 
lm 1 -0.06 


OQ + H7 
yo -76 


0.98/880 


0.85«;i 2 


l-5T±g;iS 


1 64+ - 05 


0.96/880 


PKS 0438-436 


9 10 +297 


^^-O.ll 


< 120 


0.97/280 


<0.5 




1 41+008 

r-11-0.09 


0.99/280 


PKS 0537-286 


9, +789 


1 40+ - 12 

1 -^ u -0.12 


< 4350 


0.77/203 


<1.32 


<3.8 


1 47+0.37 

-0.16 


0.81/203 


PKS 2126-158 


3981:111 


1 - oz -0.07 


< 455 


1.01/410 


1 02+ 20 


1 2fi +0 - 42 

i-^"_Q .26 


1 62+ - 11 

1 - oz -0.07 


1.00/410 



Note: Intervals are at 90% confidence, for two interesting parameters. 
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C. Complex emission models 



Object 


Bremsstrahlung model (Nn free) 
JV H kT xleJdoi 
(cm 2 ) keV 


iV H 
(cm ) 


Reflection model (Ah free) 
r R 

(=A refl /A pl ) 


Xred/dof 


S5 0014+81 


18.73+^ 


^■ L -i.§\ 


1.05/572 


27.8 -H's.r 




< 0.70 


1.04/571 


NRAO 140 


22.5±1;J 


99 Q+2-6 
z^.o_ 2 3 


0.99/881 


3i-i±!:S 


1 70+ 011 


< 0.35 


0.98/880 


PKS 0438-436 


4.63t\f 2 


r 9 7 +30 
0Z - 1 -15.7 


0.98/281 


5.86j?.£ 


1 4Q+0-26 


< 0.61 


0.97/280 


PKS 0537-286 


r\ qo + 0.83 
^•<5^_0. 60 


89.71™ 


0.80/204 


9 09+I.OO 
Z - Oz -0.72 


1 qo+0.15 

1.00_Q Qg 


< 0.35 


0.81/203 


PKS 2126-158 


6-90i j 3 7 


41.9t&° 


1.01/411 


10.412:? 


1 62+ 26 


< 0.60 


1.01/410 



Note: Intervals are at 90% confidence, for two interesting parameters. 

TABLE 7: Local Interstellar CO and IR Emission towards the Quasars 



Object 




b 


CO a 
Emission ? 


Ref. 6 


IPlOO/jm 

Emission ? 


Ref. 6 


S5 0014+81 


121.6113 


18.8020 


possible 


(1) 


No 


(6) 


PKS 0332-403 


244.7640 


-54.0749 






No 


(6) 


NRAO 140 


158.9997 


-18.7650 


Yes 


(2) 


No 


(6) 


PKS 0438-436 


248.4109 


-41.5654 






No 


(6) 


PKS 0537-286 


232.9400 


-27.2924 






No 


(6) 


S5 0836+71 


143.5408 


34.4257 


No 


(3) 


No 


(6) 


PKS 2126-158 


35.9295 


-41.8679 


No 


(4) 


Yes 


(6) 


PKS 2149-306 


17.0770 


-50.7845 


No 


(5) 


No 


(6) 



a An entry of "Yes" /"No" indicates 12 CO J=l-0 line Galactic emission has/hasn't been detected in the direction of the 
quasar (see §5. 1.1. a). 

6 References: (1) Heithausen et al. (1993), (2) Bania et al. (1991), (3) Lisz & Wilson (1993), (4) Hartmann D., private 
communication, (5) Liszt (1994), (6) Moshir (1989) 

c An entry of "Yes"/"No" indicates there is/isn't contamination by cirrus (see §5. 1.1. a). 
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Fig. 1: x 2 contour plots in the Nu-T parameter space. Full line contours represent 68%, 90% and 99% confidence 
limits obtained from ROSAT and ASCA (GIS+SIS). The 90% confidence contours obtained from fitting separately the 
GIS and SIS are indicated in dashed lines. Best-fit values are indicated by a large mark (+) for ROSAT and ASCA 
(GIS+SIS), and by a small mark (+) for GIS and SIS. The vertical line represent the Galactic absorption obtained from 
radio measurements at 21 cm with associated errors of 10 19 cm -2 if taken from Elvis et al. (1989) or with a conservative 
30% error, if taken from Dickey & Lockman (1990). 

Fig. 2: Best-fit spectra and \ 2 contour plots of the combined ROSAT and ASCA data. The data are fitted with a 
single absorbed power law model, and have been binned to a S/N ratio higher than 4 for display purposes. Note that the 
ROSAT counts (lower spectra) have been normalized by the PSPC geometric area of 1141 cm 2 . Therefore the ROSAT 
counts should be read as "normalized counts/sec/keV/cm . The contours represent the 68%, 90% and 99% confidence 
limits, and the vertical lines are, like in Fig. 1, the Galactic absorption and associated errors. For S5 0014+81, the 
energy of the Fe K edge has been indicated. 

Fig. 3: X-ray photon indices vs. (a) redshift, (b) luminosity. Stars represent the ASCA slopes and errors between 
~ 1-30 keV (quasar frame) for the 8 high-z RLQs (Table 5, fit with Nn free). Triangles represent the Ginga slopes 
and errors between ~ 2-20 keV for 18 low-z RLQs as reported by Lawson & Turner (1996). The horizontal shaded-line 
represents the approximate value < V >= 1.4 of the ~ 2-10 keV slope of the XRB (Gendreau et al. 1995). 

Fig. 4: x 2 contour plots in the A^h-F parameter space for 3C 273. Full line contours represent 68%, 90% and 99% 
confidence limits for SIS0+1. Dashed lines indicate the 90% confidence contours for SIS0 chip number 1 (SOcl) and SIS1 
chip number 3 (Slc3), fitted separately. Best-fit values are indicated with marks (+). The vertical line represent the 
Galactic absorption (full line) and associated 30% errors (dotted line) from Dickey & Lockman 1990. 
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